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We argue that specific fluctuations observed in high-energy nuclear collisions can be at- 
tributed to intrinsic fluctuations of temperature of the hadronizing system formed in such 
processes and therefore can be described by the same nonextensivity parameter q character- 
izing Tsallis statistics describing such systems (for q — > 1 one recovers the usual Boltzmann- 
Gibbs approach). It means that |<7— 1|, which is a direct measure of temperature fluctuations, 
can also be expressed by the observed mean multiplicity, {N), and by its variance, Var{N). 
This allows to deduce from the experimental data the system size dependence of parameter 
q with q ~ 1 corresponding to an infinite, thermalized source with a fixed temperature, and 
with the observed q > 1 corresponding to a finite source in which both the temperature and 
energy fiuctuate. 
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Since already some time we advocate that 
most of the single particle distributions mea- 
sured in high energy colhsions when analyzed 
by statistical hadronization model should use 
its nonextensive version (for all details concern- 
ing nonextensive statistical physics and its ap- 
plications see recent review [1], our works are 
summarized there in [2]). In such approach 
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a standard Boltzmann-Gibbs (BG) exponential 
distributions, g[E) = Cexp(— £'/T), are re- 
placed by Tsallis distribution (g-exponential) , 
hq{E) = Cq[l - {I - q)E /T]^/^^''^\ and a new 
phenomenological parameter q occurs which ac- 
counts summarily for the possible intrinsic fluc- 
tuations of the usual parameter, the temperature 
T of the hadronizing system, q—1 = Var{T) /(T)^ 
(for (? ^ 1 one recovers the usual BG approach). 
Because Var{T)/{T)'^ = 1/Cy (Cy is the heat 
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capacity under constant volume) it means that 
q — 1 = 1/Cv in fact measures the heat capac- 
ity of hadronizing system. In nuclear collisions 
in which only part of nucleus participates ac- 
tively in the collision, the above picture should be 
slightly modified by allowing for the energy trans- 
fer to / from the collision region to its surroundings 
(formed by the spectator nucleons), As result T in 
the above formulas is replaced by Tg// = To + (g— 
^)Tvisc, with To = (T) and T^sc being a new pa- 
rameter characterizing this energy transfer 0,3]. 
On the other hand, it was also shown in that 
whereas independently produced particles follow- 
ing BG distributions in energy show Poisson mul- 
tiplicity distributions, the same particles follow- 
ing Tsallis distributions instead are inevitably 
distributed according to the Negative Binomial 
formula, P{N) = T {N + k) / [T {N + l)T (k)] = 



{{N)f /{l + {N)/ky'+' with k = l/{q-l) (or 
q = l + l/k) and Var{N) = {N) + {q - 1){N)'^. 

Let us now observe that for a system with 
finite size remaining in contact with a heat bath 



\N+1 



one has, following Lindhard's approach [4], that 



Var{U) + ClVar{T) = {TfCy- 



This is a kind of uncertainty relation (in the sense 
that it expresses the truth that in the case of con- 
jugate variables one standard deviation in some 
measurement can only become small at the ex- 
pense of the increase of some other standard de- 
viation Q]). Relation ([1]) is supposed to be valid 



all the way from the canonical ensemble, 

Var{T) = and Var{U) = {TfCy, (2) 

to the microcanonical ensemble, 

Var{T) = {Tf/Cv and Var{U) = 0. (3) 

It means therefore that Eq. ([T]) expresses both 
the complementarity between the temperature 
and energy, and between the canonical and micro- 
canonical description of the system. This should 
be understood in the same way as the (improper) 
eigenstates of position and momentum appear as 
extreme cases in the quantum mechanical uncer- 
tainty relations. It is worth knowing that in 
the limiting cases of Tsallis statistics was inves- 
tigated in which q was interpreted as a measure 
thermal bath heat capacity: q = 1, i.e. canonical, 
case would correspond to an infinite bath (ther- 
malized and with fixed temperature), whereas 
q = —00, i.e. microcanonical, case would corre- 
spond to a bath with null heat capacity (isolated 
and with fixed energy). All intermediate cases 
would then correspond to the finite heat capacity 
(both temperature and energy fluctuate). 

To obtain realistic (intermediate) distribu- 
tions let us start from a system at a fixed tem- 
perature T. The standard deviation of its energy 
is 



VariU) = {Tf^ = {TfCy. (4) 



Inverting the canonical distribution gxiU) one 
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can obtain 



gu{T) = -T'—j^ gT{U')dU' (5) 

and interpret it probability distribution of 
the temperature in the system. The standard de- 
viation of this distribution then yields Var{T) = 
(T)^/Ci/. Because for a canonically distributed 
system the energy variance is VariU) = {T)'^Cv 
and because for an isolated system Var{U) = 0, 
for the intermediate case the variance (express- 
ing energy fluctuations in the system) can be as- 
sumed to be equal to 



Var{U) = {TfCv^, ^£(0,1), 



(6) 



where the parameter depends on the size of the 
hadronizing source. Inserting this in Eq. ([T]) one 
gets that q depends on ^ in the following way: 

Var{T) 1 - C 



q-1 



(7) 



Cv ■ 

It is natural to assume that the size of the ther- 
mal system produced in heavy ion collisions is 
proportional to the number of participating nu- 
cleons Np, ^ ~ Np/A. Because Cy = aNp^ we 
obtain that 



1 / Np 



(8) 



As shown in detail in [S], this is precisely what 
is observed experimentally (cf.,[T^). To recapitu- 
late what was done there: 

(i) We know that if [/ = const and T = const then 
the multiplicity distribution P{N) is Poissonian 



and we also know how fluctuations of T change 
P{N) from the Poissonian form to the NBD one. 
{ii) We want now to see how big are fluctuations 
of T in our hadronizing systems formed in a col- 
lision process. It turns out that for N = const 
we have either Eq. dH) or Eq. ([3]) depending on 
whether T = const or f7 = const. Assuming now 
validity of Eq. ([1]) not only in the above limiting 
cases but also in the general case when both the 
energy U and the temperature T fluctuate at the 
same time (fluctuations of the energy U are given 
by Eq. [6]) and corresponding fluctuations of the 
temperature T are given by Eq. ([7])). 
{Hi) Knowing how big are fluctuations of T we 
can deduce fluctuations of the multiplicity A^. 

The analysis of different a,spects of such fluc- 
tuations was performed in [3] and will not be 
repeated here. Instead we would like to bring 
ones attention to the following remarks. In the 
above formalism we have (cf. Eq. (B3) of ^ 
) that qtMU) = [/3/r(A)](/3C/)^-iexp(-/?[/) 
with fluctuations given by Var{U)/ {U)"^ = 
1/N. In fact one can now invert distribu- 
tion QT^NiU) (proceeding in analogous way as 
in Appendix B of Q]) to obtain multiplicity dis- 
tributions gT,u{^) and obtain distribution of 
temperature, gu,N{T) = -^^ Jq 9T,NiU')dU' = 
[U/r{N)]{(3U)^-^ exp{-(3U) with fluctuations 
Var{T)/{T)^ ~ Var{/3)/{(3)^ = 1/N. The fluc- 
tuations of temperature obtained this way have 
gamma distribution. Because Cy oc N one ob- 
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FIG. 1: (Color on line), (a) Prediction of Eq. (HI) ver- 
sus experimental results obtained from Var{N)/{N) 
(squares) and from the {N) (circles) - see [S!] for de- 
tails. (6) Var{N)/{N) for negatively charged parti- 
cles produced in p-l-p, semicentral C-l-C, semicentral 
Si-|-Si, and Pb+Pb collisions as function of the frac- 
tion of nucleons participating in the collision, Np/A, 
compared with our predictions (see for details). 

tains Eqs. ^ and ([3]). This illustrates that we 
can deduce fluctuations (i.e., the corresponding 
probability distributions) of any quantity out of 
(r, [/, N) provided only that the other two are 
constant. The implication of this fact are still to 
be discussed. 



Let us return again to results of [s^. The 
Eq. ([8]) can be translated to direct dependence of 
Var{N)/ {N) on the number of participants Np, 
If Teff depends on q, i.e., if T^sc 7^ and there 
is some energy transfer from the hadronization 
region to the spectators, then the above depen- 
dence is connected with the dependence of {N) 
on Np. The observed decreasing of Var{N)/ {N) 
with Np implies then the nonlinear dependence 
of (A^) on the number of participants Np. It was 
also shown in that Var{N)/{N) (xl- Np/A 
(see Fig. [TJa). In fact, also transverse momen- 
tum fluctuations defined by the measure $ be- 
have in the similar way: <I> oc 1 — Np/A. It 
means then that these fluctuations are reflecting 
not so much fluctuations of transverse momenta 
but rather fluctuations of the multiplicity. 

We would like to add here a new observation 
concerning behavior of correlations in transverse 
momenta measured in [7]. As shown in [s] one can 
connect measure of fluctuations $ with measured 
covariance in transverse momenta Cov{pTi,PTj)' 



{{N)-l+Var{N)/{N)) ■ 

■COV {pTi,PTj) ■ (9) 



Because we observe that <I> << 2y^Var (pr), one 
can write approximately that 



Cov {{pTi,PTj) ^ 
2,JVar{pT) ' 



' {N)-l+Var{N)/{N) {N)q' 

(10) 
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FIG. 2: (Color on line). Our results for correlations. 
Data are from [7]. The full line is our prediction as 
given by Eq. pT|) . See text for furtlier details. 



Using now relations derived in [3| one gets that 



Cov{{pTi,PTj] 



Np 



ill) 



This is our prediction and it is shown in Fig. [2] for 
^Np=A = 1-44 MeV and for other parameters ob- 
tained from different coniparisons with multiplic- 



ity fluctuations done in 



31], namely: (px) = 400 



MeV, b = 4.8 MeV, uq = 0.642, c = 4.1, a' = I 



and Var (px) / {pr) = 0.43. We stress here this 
fact because keeping the relevant independent 
parameters (altogether 4 of them) as free, un- 
constrained by the multiplicity fluctuations data, 
would result in a simple formula which could fit 
data more exactly. 



^yCov {pTi,PTj) 
(pt) 



Cl+C2(l-^ 



\ C3Np + Ci(l-^ 



(12) 



To close let us stress again that our approach, 
which uses the g-statistics, allows to demonstrate 
the deep connection between fluctuations and 
correlations observed in transverse momenta and 
those observed in multiplicity distributions. It 
means that they all convey essentially the same 
information on the hadronizing system produced 
in high energy heavy ion collisions. 

GW thanks the Organizers of the ISMD2009 
for their warm hospitality during the conference. 



[1] Topical Issue on Statistical Power Law Tails in 
High-Energy Phenomena, Ed. Biro T S (2009) 
Eur. Phys. J. A 40(3). 

[2] Wilk G, and Wlodarczyk Z (2009) Eur. Phys. J. 
A 40: 299 - 312. 

[3] Wilk G, Wlodarczyk Z (2009) Phys. Rev C 79: 
054903(10pp) 

[4] Lindhard J, 'Complementarity' between energy 
and temperature, in The Lesson of Quantum The- 
ory, Eds. de Boer J, Dal E, Ulfbeck O (1986) 



(North- Holland, Amsterdam). 
[5] Uffink J, and van Lith J (1999) Foundations of 

Phys. 29: 655 and Thermodynamic uncertainity 

relations, cond-mat /9806 1 02[ 
[6] Campisi M (2007) Phys. Lett. A 366: 335. 
[7] Adams J et al. (STAR CoUab.) (2005) Phys. Rev. 

C 72: 044902 (6pp). 
[8] Utyuzh O V, Wilk G, Wlodarczyk Z (2001) Phys. 

Rev. C 64: 027901 (4pp). 



5 



